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I. INTRODUCTION 



Recently, substantial progress has been made in developing Chiral perturbation theory 
(ChPT) for pionic nucleon (irN) system □ □ □ D □ □ □ In this context the precision 
X-ray experiment on pionic hydrogen and pionic deuterium carried out by ETH Zurich- 
Neuchatel-PSI collaboration 8}. The s— wave scattering nN scattering lengths are important 
for testing of various theoretical consideration like the Goldberger-Miyazawa-Oehme sum 
rule |9j determining the ttttN coupling constant. The sigma term, which is used for lattice 
and ChPT calculation, is sensitively affected by the isoscalar scattering length. The tcN and 
7i d scattering lengths are determined as directly from the phase shift analysis so from the 
X-ray experiments on pionic atoms using the Deser formula lo| . 

The extraction of the s-wave scattering length from the X-ray experiments gives errors 
less than the phase shift analysis. The accuracy of the modern level of experimental analysis 
of the parameters of pionic hydrogen reached by the PSI Collaboration is about 0.2% for 
the strong energy shift and 1% for the width of the energy level of the ground state of pionic 



hydrogen 11 1. The NPE measurement |l2| of the pionic deuterium atomic level strong shift 
yields for the ird scattering length a nd = [-0.0261(±0.0005) + «0.0063(±0.0007)]m- 1 , where 
m-n is pion mass. Note an imaginary part of ltd scattering length approximately in four times 
less than its real part. Such results for ix~d scattering can be explained by fact that the 
absorption channel n~d — > nn and the radiative absorption channel ir~d — *> ^nn give very 
small contribution to the scattering length. Also ir~d — > 7r°nn channel, which is opened at 
threshold, is suppressed by the centrifugal barrier. Indeed, the 96% of the deuteron state is 
3 S\ state, therefore according to Pauli principle 3 Si is forbidden for nn system and 1 S is not 
available if there is no spin flip. In the series of papers 13|, llJ, Ilia, [16J the determination of 



the pion-deuteron scattering length was considered by solving the Faddeev equations. There 
has been remarkable recent progress in developing the effective field theories to problems 
relevant to the nN jl^, ^| and 7id j^J scattering length. The main results which obtained 
from calculation of the nd scattering length are following: (a) in all theoretical calculations 
isospin symmetry of the strong interactions is assumed; (b) the ir~d scattering length can 
be extracted from the Is energy shift of pionic deuterium using the Deser formula; (c) the 
electromagnetic contribution to the difference of the real parts of the scattering length of 
7r + and 7r~ is of order of the experimental error on the scattering length. The aim of this 
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paper is to establish the precise relation between the strong energy-level shift of the 7r~d 
atom in the Is state, and strong ird scattering length using the effective field theory. In the 
present paper we show how to obtain the nonrelativistic Hamiltonian for ird system from 
the relativistic equations of motion for the scalar and vector fields with taking into account 
the electromagnetic interaction. The strong interaction is taken at the zero-range approach. 
Here we use the unit system h = c = 1, therefore the fine-structure constant is a = e 2 . 



II. NONRELATIVISTIC LIMITS FOR THE SCALAR AND VECTOR FIELDS 

The problem of nonrelativistic limit description for fundamental particles and their inter- 
actions may be solved in different ways. Although in all methods of nonrelativistic expansion 
the first terms of the Hamiltonians coincide, however the difference begins to arise at tran- 
sition to the higher orders of expansion. The method of Foldy-Wouthuysen transformation 
is one of the safest method for solving this problem j^jj. This method was mainly used for 
spinor and scalar particles. In our task we need in nonrelativistic Hamiltonian for massive 
vector field (deuteron is considered as a fundamental particle) interacting with electromag- 
netic and scalar (the pion) fields. In this approach suggested by us the equations of motion 
are basic for producing nonrelativistic Hamiltonian. This Hamiltonian we obtain in the form 
of series over 1/M, where M is the mass of our vector particle. 

We would like to point out some steps of derivation of the nonrelativistic limit Hamilto- 
nians by means of the Foldy-Wouthuysen transformation [2l[. The method is based on the 
transformation of a relativistic equation of motion to the Schrodinger equation form. For 
example, the relativistic equation of motion for free massive scalar field ip is 

(d 2 + m 2 ) V = 0, (1) 



or 



d\p> = (V 2 - m 2 ) <p, (2) 

where m is mass of the scalar particle, d 2 is the four-dimensional D'Alambertian and dt = 
d/dt (we use metric goo = +1 , c?ii = 922 = #33 = —1)- Let's introduce the "big" 6 and the 
"small" x components of the field ip as 

1 i 1 i 

e = ^ + — dtp), X = t;(<P %>)• ( 3 ) 

2 m 2 m 
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In the case ip ~ e imt (i.e., when ip corresponds to a particle) we have 9 ~ ip, x ~ 0, and 
vice versa in the case <p ~ e +* m * W e have 6 1 ~ 0, x ~ V 9 - ^ is easy to see that 



V 2 

id t 8 = m9-— (6 + x) 
2m 



V 2 

idtX = ~ m X + 7^—( 6 + XJ- 
2m 



(4) 



If we introduce a new two component function as 

-(!)■ (5) 

then we can rewrite Eq. (J2J) in the Schrodinder equation type 

id t $ = H®, (6) 

where 

E = J m + (7) 



and 





fl ^ 




r l\ 


v = 1 




1. H 


I -1 J 




10-1/ 







(8) 



The matrix p is that one which mixes the "big" and "small" ("particle" and " antiparticle" ) 
components of $. Using the unitary transformation 

$' = e l5 $, #' = e iS He~ iS , (9) 

searching the operator 5 in the form 

S = ijp\, (10) 

and demanding removing of terms of odd order in p in the Hamiltonian we find the operator 
A of the Foldy-Wouthuysen transformation to be 

*tan(2A) = ^P^, (11) 
m l + Pq 

and the new Hamiltonian 



H' = t/Pq = riym 2 + p 2 . (12) 
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So, our new Schrodinger equation idt& = H'<f>' is free of mixing terms. 
Now we include the electromagnetic field. The equation of motion is 

(D 2 + m 2 )ifi = 0, or D%(p = (D 2 - m 2 ) ip, (13) 

where we have used the standard notations 

Dp — dp — ieA^, D — d — ieA , D = -V - zeA, 7r = iD = p + eA, (14) 

and is the electromagnetic potential, 7r is generalized momentum. Determining the "big" 
and "small" components of (p by equations 

6> = v?H D <p, x = <-P D tp, (15) 

m m 



we receive 



or 



D 2 D 2 

tD 6 = m8- — (8 + X ), iD oX = -m X + ^( 6 + x), (16) 



where 



id t $ = H<f>, (17) 

D 2 D 2 7T 2 7T 2 

H = mt] - eA - 77- p— = mr] - eA + r]- hp—. (18) 

2m 2m 2m 2m 

This Hamiltonian we should transform to separate the states with the positive energy from 
that of negative ones. Substituting 



$ = e - iS & (19) 



into Eq. (fTTj) we obtain 



where 



i3t$' = (20) 



H' = e lS (H - d t S) e~ iS . (21) 
Using standard expansion formula 



e iA Be -iA = B + t[AjB] + ^ [A} [Ai b]} + t[A, [A, [A, B}]] + ... (22) 



and taking into account that S ~ 1/m, which follows from the free case, we get 

H> = H + i[S, H] - hs, [S, H]] - *-[S, [S, [S, H]]] 

2 o 

+ ±[S, [S, [S, [S, H]]]] - d t S - d t S] + ±[S, [S, d t S\], (23) 



up to the accuracy 0(l/m 4 ) [21 1. 



Let us use the following notation 



where 



D 2 D 2 

H = mil — eA — r> p — = rim + £ + 1Z, (24) 

2m 2m 



7T 2 7T 2 

S = -eA Q - i]B 2 /2m = -eA + n — , K = -pD 2 /2m = p — . (25) 

2m 2m 



We suppose that the Coulomb energy has the same order of magnitude as the kinetic energy. 
It is usual situation for nonrelativistic case. Our aim is to eliminate terms with odd degrees 
of 1Z. This work may be done only step-by-step. After some steps we obtain the expression 
which is free of odd degrees of 1Z with sufficient accuracy (see |2l|): 

H' = V (m+^-TZ 2 - -^K 4 ) +£- [K,£]) - -^[K, d t K] + 0(l/m% (26) 

\ 2m 8m d J 8m 2 8m/ 



It is not difficult to show that 



Z6 

[R, £] + id t K = P 7v ■ E. (27) 

m 



So, we have the transformed Hamiltonian up to order 0(l/m 4 ) in the form 



2 4 , 

7T 7T \ , ie 



^^ '=^ ro+ ^"8^J" cAo "^ l "'' r ' E, ' (28) 

Here we have no terms mixing the "small" and the "big" components. This Hamiltonian 
includes terms of order (v/c) 2 . 

Now we can proceed to a massive vector field theory. A massive particle of spin 1 has 
hxee degrees of freedom, so , it may be considered as spatial part of a four vector tp^. In Ref. 



22| a general theory describing particles of unit spin and arbitrary magnetic moment was 



developed and applied to the motion of such particles in the electromagnetic field. For our 



24j . The massive four-dimensional 



purpose we use the equations of motion given in Ref. [23, 
vector field ip u is described by equations 

VV = D^u - D v %- + - ie/cF^ = 0, (29) 

where M is mass of the vector particle, 

Fun = d ^ - d »K (3°) 
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is the strength tensor of the electromagnetic field, 

D„ = - ieA^ (31) 

is the electromagnetic covariant derivative and k is the anomaly part of magnetic moment 
of the particle. Acting on the second equation by D v we receive the following constraint: 

^„ = -^^vv + ^(F^r). (32) 

This constraint reflects the fact that in the Lorentz covariant description of a massive vector 
field with three degrees of freedom we have used a four vector ip^, so one of its components 
must be excluded. 

Let us rewrite the Eq.([29j) as 

D 2 ip u - D^D^ + M 2 4> u - ienF u ^ = 0, (33) 



D 2 ipv ~ D U D^ + M 2 i) v - ze(l + k)F v ^ = 0, (34) 



or, 



Using the constraint (|32|) we get 

Dfy v = (D 2 - M 2 ) i, v + ie(l + k)F v ^ - — [D^F^r,) - 2kD u D»(F^ x )] . (35) 

When it needed we should replace the zero component of ip v through its spatial part (cf. 
Eq. (HD ): 

V>o = M2 ^ p2 [ie(l + «)E ■ V + A)£>VJ , (36) 

where ipi are spatial components of if). It is not difficult to convince that Eqs. ()32j) and (jHEj) 
coincide. Further we will write out equations for the spatial part of only. 
Let us to introduce the following notations 

% % 



Then 



iDati = iDoik - ^D&t = MQ - - + (38) 

iD m = iD^i + ^D 2 ^ = -M Xl + ^ + W{1 ^ K) F lX ^ x - O(j^). (39) 
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Supposing 



Xi 



we can rewrite the last system of equations as 

"VI', (41) 
Nonrelativistic equations for components of may be written as 

id Q — (M — eAo)d - ^(Q + Xi) + (B x (C + X )\ + 0(±) , (42) 

id oXi = — (M + eA ) Xl + ^(d + X.) - (B x (C + x)h + 0(±) ■ (43) 

In the two-component notations we have 

= (-^o + M V - + p)) % + ^^-(v + P)(B x *), + 0(±). (44) 

The situation is like to the scalar case except the term with magnetic field. Denoting 

Hij = Mr]5 tj + Sij + Kij (45) 



we obtain 



(D 2 \ ie(l + k) 

~ CA ° ~~ V 2M J 5ij + 2M v£ik i Bk > ( 46 ) 
m D 2 p ie(l + k) 

ij = ~ P 2M ij + 2M p£ik J Bk - ( A7 > 



Let's at the first step restrict by terms of 1/M order. In this case we have 



(48) 



For the Hamiltonian density we receive 



v { m+ ^m)- eA °\ + ^W^* dV Bx * d > (49) 

where Sf^ denotes the non-relativistic vector field operator, are components of This 
is a nonrelativistic Hamiltonian where the spin (magnetic moment) of ^-field is taken into 
account. 

The procedure of getting the next terms is the same as for scalar field. The result of the 
Hamiltonian without the terms mixing particle-antiparticle is: 

= " ( M + 535?** - sW nl ) + £ - sW^' ^ f « - slr-^ W < 50 > 



Here it is assumed, that all terms enter with the additional tensor structures which are 
defined in Eqs. (gSJ) - (EE7I) 

The terms up to order 1/M 4 are included into Eq. (J5U)) . In this work we will restrict 
ourself by terms of order 1/M 3 only omitting all terms of higher accuracy. Therefore for the 
vector field we obtain the Hamiltonian density 



■q [M + 



7T 



P 1 



2M 8M 3 



eAr 



** + B >< 



(51) 



Now we can write the total Hamiltonian including into consideration other fields that 
means pion and electromagnetic fields. Also we include into the density of the Hamiltonian 
the strong 7rd interaction which is taken in the zero range approach 



(52) 



where <&„- denotes the non-relativistic field operator of pion and d^d is the coupling constant 
which does not depend on the spin projection of deuteron. The electromagnetic Lagrangian 
is 



_V „ = -(E 2 -B 2 ) 



-E • «9 A + AjdivE (E 2 + B 2 



i.e., 



1 



(53) 



(54) 



is the (density) of the electromagnetic Hamiltonian and A Q is a Lagrange multiplier (it is 
not a dynamical variable). 

Taking into account the Lagrange multiplier term we get the total Hamiltonian density 
in form 



TjlM + 



7T~ 



2M 8M 3 



+ eA c 



+ -(E 2 + B 2 )-A)divE + $t 



7] I m + 



di 



7T" 



ie(l + k) t + n T 



2M 

,-.4 



2m 8m 3 



(55) 



In the gauge divA = we have divE = — AA (due to E = — VA — d t A), so, the Coulombic 
part of our Hamiltonian is 

A AA + ^(VA ) 2 + pA , p = e(^ d -^ d - 
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Excluding Aq and neglecting the spin-flip part of the density Hamiltonian (|55j) we obtain 
the canonical Hamiltonian for the non-relativistic nd system as 



H = H + H c + H R + H 7 + H s = H + H c + V, (56) 

where Ho is the free Hamiltonian describing non-relativistic pions and deuterons. Further, 
H r = / d 3 xW r , T = C, R,7, S, and 

H c = e 2 (*;*,) A" 1 (d^) , (57) 

? f IP 

H ~< = ~2M * d (VA + AV) *" + 2^ (VA + AV) ^ ' (59) 

H s = -d^Vfod^. (60) 



III. ENERGY SHIFT OF PIONIC DEUTERIUM ATOM 

We consider the problem of s-state energy shift according to the perturbation theory. 
Such analysis was performed for the pionic hydrogen in Ref. [18] . Let Ho + He is the 
unperturbed Hamiltonian, whereas V is considered as a perturbation. The ground-state 
solution of the unperturbed Schrodinger equation in the center of mass (CM) system frame 
(Ex - H - Hc)|^o(0)) = 0, with Ex = M + m + Ex, is given by 

l*o(0)> = / S3 *o(p) ^(p) at(-p) |0> , (61) 

where a'(p) and b'(p) denote creation operators for non-relativistic tt~ and deuteron acting 
on the Fock space vacuum, and ^o(p) stands for the non-relativistic Coulomb wave function 
of IS state in the momentum space, E\ is the non-relativistic binding energy of pionic 
deuterium. We normalize the ground-state wave function as 

J d 3 p|^ (p)| 2 = l. (62) 

We are going to evaluate the energy-level shift of the ground state due to the perturbation 
Hamiltonian V as in Ref. Let us define the free and Coulomb Green operators by the 
expressions G (z) = [z — H ) _1 and G(z) = (z — H — H c ) _1 , respectively. Further, 
we define the "Coulomb-pole removed" Green function as G(z) = G(z)(l — U), where n 
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denotes the projector onto the Coulomb ground state \l/o (EU)- The ird scattering states in 
the sector with the total charge are defined as |P, p) = ^(pi) <^(p2) |0) (pi and P2 denote 
momenta of deuteron and pion, respectively). The CM and relative momenta are defined by 
P = pi + p 2 , p = {mpi — Mp 2 )/(M + m). We remove the CM momenta from the matrix 
elements of any operator R(z) by introducing the notation 

(P,q|R(z)|0,p) = (27r) 3 5 3 (P)(q|r(^)|p). (63) 

The "Coulomb-pole removed" transition operator satisfies the equation 

M(z) = V + VG(z)M(z) . (64) 

According to the Feshbach's formalism j2] the scattering operator T(z) develops the 
pole at z = z where z is the solution of the following equation 

f-^ 1 -(*o|m(f)|* o ) = 0, (65) 

where (p|\l/o) = ^o(p) an d m(z) is related to M(z) through the definitions (|63|). 

In order to get the shift of the ground-state energy, the quantity m(z) is calculated 
perturbatively from Eq. ()64)1 by the iteration series at accuracy 0(a 4 ). Using the explicit 
expression of V given by Eq. (J56|) . replacing G(z) by Gq(z) whenever possible, and retaining 
only those terms that contribute at the accuracy we are working, the operator M(z) can be 
written in the form M(z) = U(z) + W(z), where 

U(z) = H R + H 7 G (z)H 7 , 

W(z) = H s + H s G(z)H s , (66) 

At the accuracy 0(a 4 ), the energy of the bound state is equal to z = E\ + AE° m + ei s , 
where 

AEr = Re(tf |u(£i)|tfo) + } e ls = Re (tt |w(£i)|tt ) ■ (67) 
Here u(z), w(z) are related to U(z), ~W(z) through the definitions ()63|) and E™ c stands for 
the contribution due to the electron vacuum polarization which is added "by hand". The 
results of calculations for the matrix elements determining the energy-label shift are the 
same as in Ref. jl^ : 

R.(».Wft)W = -^^-^, (68) 



Re(#o|w(Ei)|tf ) = 



7T 



2 

_d jrd + ^(e+^(ln Qt -l)) 

7T 
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Table 1. Contributions to the electromagnetic binding energy of the it p and it d atoms (eV). 



Type of contribution Notation tt p tt d 

Point Coulomb, KG equation E KG -3235.156 -3459.0 

Vacuum polarization, order a 2 E vac —3.241 —3.732 

Relativistic recoil E rel 0.034 0.021 



where fj, % a is the reduced mass of the ird system, and 



a f J 'TTd I / 2\d— 3 / -p'/i \ 1 A I i 1 



£ = ^^(^r 3 (T^- r, (l)-ln47r)+ln^^-U. (69) 



2tt [ xr ~ ' \d-3 ~ v ^ J ' ~ /i 2 

Here d and // denote the dimension of space and the scale of the dimensional regularization 



used as in Ref. [18|], respectively. 

The energy shift ( order a 2 ) due to the vacuum polarization contribution is given by the 
well-known expression jl^ 

£ vac = -^V$(??), U = ^, (70) 

r? 2 $ (77) = - 2tt - 4r? + — r? 2 - — r/ 3 + ' In (rj + vV - 1) , 

T\ 2 3 J rf^Ji] 2 - 1 

where m e denotes the electron mass. 

The calculation of the electromagnetic energy-level shift is now complete. 
We present our results in the convenient form 

E x + AEl m = E KG + £ vac + E rel , (71) 

where 

E KG = + (72) 

E rel = l a ±A± (73) 

8 Mm ' v ; 

and E vac is given by Eq. (|70p. In Table 1 we present various contributions to the electromag- 
netic energy-level shift of 7i~p and ir~d atoms. As can be seen from Table 1, the calculated 
values of the electromagnetic shifts for ir~p and 7T~d atoms are close to the each other. 
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IV. S- WAVE SCATTERING LENGTH AND THE DESER'S TYPE FORMULA 



In order to complete the calculation of the strong energy- level shift, one has to match at 
the accuracy 0(a) the particular combination of the non-relativistic coupling d^d- Therefore 
we consider the scattering operator 

T r (*)=Vr + V r G r (z)T r (z), 

V R = H c + H 7 + H s , G R (*) = (z - H - Hr)" 1 . (74) 

In the scattering operator Tr(z), all kinematical insertions contained in Hr are summed 
up in the external lines (see U for the details). We calculate the matrix element of the scat- 
tering operator Tr(z) between the 7r~ei states at 0(a). After removing the CM momentum, 
the spin-nonflip part of this matrix element on energy shell is equal to 

(qltaMlp) = - + e**<« (qfeWlp) , 

|q — p\ 2 4Mm |q — p\ 2 

(75) 

where the (divergent) Coulomb phase in the dimensional regularization scheme is given by 

fc(|p|) = ^f/- 3 ( ~ + 1*4,0) + In M , (76) 
|p| yd — 3 2 / |p| n 

and 

Re (q|t a (s)|p) = _Z^k + ln IPI _ ^ + dU + . . . ( (7?) 

|P| 7T /i^ 

where ellipses stand for the terms that vanish at threshold. The first two terms in Eq. (|77|) 
are arisen due to the electromagnetic interactions between deuteron and pion, the sum of 
two next terms is amplitude of scattering due to the strong interaction at threshold. 
Therefore for the regular part of the s-wave scattering length we have 

2tt 

A^d = -dvd + dl£. (78) 

Finally compare Eq. ()78|) with Eq. where we take into account all terms at accuracy 
0(a 3 ) we obtain the relation between the regular part of the s-wave nd scattering length 
and the strong energy- level shift of n~d atom: 

ei, = -2aV^^rd, (79) 
where the ultraviolet divergence contained in the quantity £, has been canceled. 
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V. DISCUSSION AND CONCLUSION 



From the measurement of the 3p — Is X-ray transition of pionic deuterium the n d 
scattering length has been determined. The strong energy-level shift can be defined by 

e ls = Ega u - E™1 , (80) 

where is the electromagnetic transition energy calculated in the absence of the strong 

interaction [27[ (the strong interaction in the 2>p state is negligible), E^ e _^( s is the measured 
transition energy. It is clear that on the extracted significance of €\ s energy shift the con- 
tributions to the electromagnetic energy E^°f ls are affected. According to our suggested 
model the contribution of deuteron magnetic moment would be differed from the results of 
the standard model using for calculation of the electromagnetic transition. Therefore it will 
be useful to estimate this contribution. 

In order to combine the results from pionic hydrogen and pionic deuterium the establish 
of relation between A n d and ttN isoscalar and isovector lengths &o and b\ is the important 
problem. This relation can be expressed as 

1 +m/m N (higher order) ( * 

= YT^/W nd ' (81) 

where the first terms is the scattering amplitude in the impulse approximation, 26o is the 
sum of the amplitude for n~p and n^n elastic scattering, is nucleon mass. The second 
term is dominant and it includes all remaining (higher order) contributions. The results 
presented in Ref. js| shows that the s-wave double scattering gives the most contribution 
while the multiple scattering and form factor correction are much smaller. It is well-known 
that the second-order formula for ird scattering length is described by (cf. Ref. [28]) 



in 



bo + (K - K) 



(82) 



where the expectation value of 1/r is taken with respect to the deuteron wave function and 
bj = (1 + m/mAr)6o- The calculations of the expectation value performed by the various 
NN potentials gave 20% of discrepancy for its value Q]. That means it is necessary to 
calculate the expectation value with potential which has the right behavior at the small 
distance. We note that absolute value of A n d (Eq. (fS2|0. calculated with use of various NN 
potentials leads to the larger value than at the solution of the Faddeev equation with the 
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same potentials. On the basis of the solution the Faddeev equation received in Ref. 
formula (jH2j) can be extended 

-2nr 



the 



•And 



m 



bo + (K - K) 



+ (bo + b 1 ) 2 (bo-2b 1 ] 



(83) 



where n = ^2fi np €d and is the binding energy of the deuteron, is reduced mass for 
proton and neutron. In contrast to Eq. (}%2*j) the binding energy correction of deuteron and 
the other higher order corrections are taken into account in Eq. (jH^j) . 

At the end we would like to emphasize while we have not taken into account the im- 
portant parameters of deuteron andpion as their sizes, but it is possible to construct the 
non-relativistic Lagrangian, like in 18|, which includes these parameters through coupling 
constants. 

We are grateful to Prof. J. Gasser and Dr. A. Rusetsky for the current interest in 
the work and useful suggestions. This work was supported by the Swiss National Science 
Foundation, SCOPES 20002003, Project No. 7UZPJ65677. 
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